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Electric charge fluctuations 
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Results are presented on event-by- event electric charge fluctuations in central Pb+Pb col- 
lisions at 20, 30, 40, 80 and 158 AGeV. The observed fluctuations are close to those 
expected for a gas of pions correlated by global charge conservation only. These fluctua- 

>^ ! tions are considerably larger than those calculated for an ideal gas of deconfined quarks and 

gluons. The present measurements do not necessarily exclude reduced fluctuations from 

2 ■ a quark-gluon plasma because these might be masked by contributions from resonance de- 



cays. 
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I. INTRODUCTION 

Ultra-relativistic heavy-ion collisions provide the opportunity to study the properties of 
strongly interacting matter. One of the predicted features of this matter, which one hopes 
to establish in heavy-ion collisions, is the occurrence of a phase transition between a purely 
hadronic state and the quark-gluon plasma. Recently several results were reported [l|, l^l 
which suggest that this transition starts in central Pb+Pb collisions at energies around 30 
AGeV J2I y . The search for further signals of deconfinement is in progress and may provide 
additional support for such an interpretation. Among them a suppression of event-by- 
event fluctuations of electric charge was predicted J5|,|6| as a consequence of deconfinement. 
Estimates of the magnitude of the charge fluctuations indicate that they are much smaller 
in a quark-gluon plasma than in a hadron gas. Thus, naively, a decrease of the fluctuations 
is expected when the collision energy crosses the threshold for the deconfinement phase 
transition. However, this prediction is derived under the assumptions that the initial 
fluctuations survive hadronization and that their relaxation times in hadronic matter are 



B 



significantly longe. than the hadronic stage of the co,U.on BHH. F« data on change 

fluctuations in central heavy ion collisions were published by PHENIX p] and STAR 
at the BNL RHIC, and preliminary results at the CERN SPS were shown by NA49 [lQ| . 
The predicted large suppression of charge fluctuations was not observed. Results by N A49 
on transverse momentum and strangeness fluctuations can be found in Refs. [ll|, |l2| and 



m 



IJ], respectively. 

In this work final results on electric charge fluctuations in central Pb+Pb collisions at 
20, 30, 40, 80 and 158 ^GeV measured by NA49 at the CERN SPS are presented and 
discussed in view of their significance as a signal of deconfinement. The used measure of 
charge fluctuations A<l>q |l^ is introduced in Sec. HH The experimental set-up is presented 
in Sec. Illll and the data sets as well as analysis method are described in Sec. IIVI Results 
are given in Sec. V and are discussed in Sec. VI. The summary is given in Sec. VII. 

II. THE MEASURE OF CHARGE FLUCTUATIONS 

The magnitude of the measured charge fluctuations depends not only on the unit of 
electric charge carried by degrees of freedom of the system (hadrons or quarks and gluons) , 



but depends also on trivial effects, which may obscure the physics of interest. The two 
most important of these effects are the fluctuations in the event multiplicity, caused mostly 
by the variation of the impact parameter, and changes in the mean multiplicity due to 
changes of the acceptance in which fluctuations are studied. In addition to the D measure 
of cKa.,e fluc.uaUo. Q sev.a. a.en.«ve .eas.e. such as .,.,,„, Q a.d A*, Q 
were proposed to minimize the sensitivity to these effects. In this analysis we use A<l>g 
which is constructed from the well established measure $ of event-by-event fluctuations, 
defined as |l8[: 



where: 



TV 

z = X — X, Z = y^{xi — x). (2) 

In these equations x denotes a single particle variable, A^ is the number of particles of 
the event within the acceptance, and over-line and (...) denote averaging over a single 
particle inclusive distribution and over events, respectively. By construction, for a system 
which is an independent sum of identical sources of particles the value of $ is equal to 
the value of $ for a single source and does not depend on the number of superimposed 
sources (l8l.ll9l|. In this analysis x in Eqs. (jS)) is taken to be the electric charge q and the 
measure is called $g. 

For a scenario in which particles are correlated only by global charge conservation 



(GCC) the value of $g is given by 
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*,,Gcc = Vr^i^-l, (3) 

where 

P=J^^ (4) 

{Nch)tot 

with {Nch) and {Nch)tot being the mean charged multiplicity in the detector acceptance 
and in full phase space (excluding spectator nucleons), respectively. Strictly speaking 



Eq. (3) holds for vanishing net charge. However, as shown in 



3, 



Eq. (3) serves as a 



good approximation also for realistic non-zero values of the total net charge. 



In order to remove the sensitivity to GCC the measure A$„ is defined as the difference: 



A$. = $. 



$n 



(5) 



'q — ^q ^q,GCC ■ 

By construction, the value of A$g is zero if the particles are correlated by global charge 
conservation only. It is negative in case of an additional correlation between positively 
and negatively charged particles, and it is positive if the positive and negative particles 
are anti-correlated Il5ll. 



III. EXPERIMENTAL SET-UP 
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FIG. 1: The experimental set-up of the NA49 experiment. 



The NA49 experimental set-up J20| is shown in Fig. Q] The main detectors of the 
experiment are four large- volume Time Projection Chambers (TPCs). Two of these, the 
Vertex TPCs (VTPC-1 and VTPC-2), are located in the magnetic field of two super- 
conducting dipole magnets. This allows separation of positively and negatively charged 
tracks and a measurement of the particle momenta. The nominal magnetic field is adjusted 
in proportion to the beam energy to ensure a good acceptance at all energies. The 
other two TPCs (MTPC-L and MTPC-R), positioned downstream of the magnets, are 
optimized for precise measurement of the ionization energy loss dE/dx which is used for 
the determination of the particle masses. Additional information on the particle masses 
is provided by two Time-of-Flight (TOF) detector arrays which are placed behind the 
MTPCs. The centrality of the collisions is determined by a calorimeter (VCAL) which 
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measures the energy of the projectile spectators. To cover only the spectator region the 
geometrical acceptance of the VCAL was adjusted for each beam energy by a proper 



setting of a collimator (COLL) |20|, |2l|. The Beam Position Detectors (BPD-1, BPD-2 
and BPD-3) are used to determine the x- and y-coordinate of the beam at the target. 
Alternatively the main vertex position is reconstructed as the common intersection point 
of reconstructed tracks. A detailed description of the NA49 set-up and tracking software 
can be found in Ref. 

IV. DATA ANALYSIS 



Data 



At all energies 50K events were analyzed with a centrality of 7% of the inelastic cross- 
section except at 158 AGeV where the 10% most central events were selected. 
To minimize the contributions of non-target collisions only events which satisfy the fol- 
lowing two selection criteria were used in the analysis. Firstly, the reconstruction of the 
primary vertex position based on BPD and TPC data had to be successful. Secondly, the 
difference between vertex coordinates resulting from the BPD and TPC data should be 
smaller than ±1 mm in x- and y-coordinate and ±5 mm in ^-coordinate. 

Several quality criteria were applied to the particle tracks. All tracks should contain 
points measured in at least one of the Vertex TPCs and the number of all measured 
points, np, should be larger than 30. Only for these particles charge and momentum 
determination is considered to be reliable. To avoid double counting of particles only 
tracks with a measured number of points larger than 50% of all geometrically possible 
points were accepted. The number of particles originating from weak decays and secondary 
interactions is reduced by only using tracks for which the x- and y-position extrapolated 
to the z-coordinate of the target is close to the position of the interaction point (Ife^^l < 2 
cm for the x- and \by\ < 1 cm for the y-coordinate) . 

Furthermore particles are required to lie in a well defined acceptance region in y, pt and 
(f) {y is the rapidity in the center-of-mass system calculated assuming the pion mass, pt is 
the transverse momentum and denotes the azimuthal angle). A well defined acceptance 
is essential for comparison of the results with model predictions and with data from other 



experiments. The acceptance limits are parametrized by the function: 

1 



PT{y,(f)) 



My) + {^r 



+ B{y), 



(6) 



where A{y), B{y), C{y) and D{y) are parameters depending on the rapidity and colhsion 
energy. The values of the parameters for positively charged tracks in the nominal magnetic 
field {By pointing upward) are given in Table 1. These parameters also apply to negative 
tracks or to a reversed magnetic field {By pointing downward) provided in Eq. (6) is 
replaced by cp' = sign{(f)){180 — \(f)\). 

As an example we show in Fig. |21 the acceptance in p^, used for the analysis for 
-0.2 < y < and 1.4 < y < 1.6 at 20 ^GeV and -0.6 <y < -0.4 and lA < y < 1.6 at 
158 v4GeV. 
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FIG. 2: Distributions of registered particles in the pr-^-plane for —0.2 < y < (a) and 
1.4 < y < 1.6 (b) at 20 ^GeV and -0.6 <y < -0.4 (c) and 1.4 < y < 1.6 (d) at 158 AGeV. 
The acceptance limits used in the analysis are shown by the solid lines. 



TABLE I: Values of the parameters A, B, C and D of the acceptance Hmits, Eq. 6, for different 
energies and rapidities. The dimensions of the parameters are such that the use of the (p angle 
in degrees and the rapidity in the center-of mass system results in the px limit in GeV/c. 
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B. Analysis 



Charge fluctuations are studied as a function of the width of the rapidity interval Ay. 
These rapidity intervals were centered around 1.27, 1.07, 0.89, 0.89 and 0.89 for the 20, 



30, 40, 80 and 158 AGeV data, respectively. The measures $g and A$g were calculated 
for ten different rapidity intervals increasing in size from Ay = 0.3 to Ay = 3 in equal 
steps and will be plotted either versus Ay or the corresponding ratio {Nch)/{Nch)tot The 
largest rapidity interval contains approximately 90% of all accepted particles. In Fig. 3 
the rapidity distribution of the accepted particles at 158 AGeV is shown together with 
the largest rapidity interval used in the analysis. 
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FIG. 3: The rapidity distribution of accepted particles in central Pb+Pb collisions at 158 ^GeV. 
The largest rapidity interval used for the analysis is indicated by dashed lines. 



For each event the positively and negatively charged particles which fall into each rapidity 
interval and the corresponding Pt-4> acceptance are counted and using these numbers 
(A^+ and A^_) the values of A$q are calculated. The total charged particle multiplicity, 
{Nch)tot, was estimated for each energy based on the NA49 measurements |l|, |^. 
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Errors 



The statistical error of A$q is calculated by dividing the whole sample of events into 
10 subsamples and calculating A$g for each subsample separately. The dispersion of 
the obtained A$g values divided by -\/9 has been taken as the statistical error. The 
systematic errors of A$g are estimated by varying track quality cuts: The values of A$g 
are calculated for two additional sets of cuts, more {np = 35, \bx\ < 0.75 cm and \by\ < 0.5 
cm) and less {np = 30, |6^| < 4.5 cm and \by\ < 2.5 cm) restrictive in comparison to the 
standard cuts. The accepted particle multiplicity decreases by about 25% when changing 
from less to more restrictive cuts. The difference of these two A$g values is considered 
as the systematic error. Since the statistical errors are much smaller only the systematic 
errors are shown in the figures. 

V. RESULTS 

A simple measure of charge fluctuations is the width of the distribution of net-charge 
Q = N^—N- in the events. As an example the distribution for central Pb+Pb collisions at 
158 AGeV is shown in Fig. |3] This distribution is compared to the net-charge distribution 
obtained from mixed events (dashed line in Fig. |3]) constructed by randomly selecting 
particles from different events according to the multiplicity distribution measured for the 
data. The net-charge distribution from mixed events is significantly broader than the 
net-charge distribution obtained from real events. 

The main source of the observed difference is charge conservation which correlates 
positively and negatively charged particles in the real, but not in the mixed events. This 
is demonstrated in Fig. El where the $g values are plotted as a function of the fraction of 
accepted particles {Nch)/{Nch)tot for central Pb+Pb collisions at 20-158 ^GeV. The main 
trend observed in the data, a monotonic decrease with increasing fraction of accepted 
particles, is approximately reproduced by introducing global charge conservation as the 
only source of particle correlations (Eq. 3 shown by the dashed line in Fig. Ej). 

By construction, the previously introduced measure A$^ is insensitive to the correla- 
tions due to global charge conservation (see Sec. HH). The dependence of A<l>g on the width 
of the rapidity interval Ay is shown in Fig. El for central Pb+Pb collisions at 20, 30, 40, 80 
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FIG. 4: The distribution of the net-charge for central Pb+Pb colUsions at 158 ^GeV (sohd hne) 
and the corresponding distribution obtained for mixed events (dotted Une) in the maximum 
rapidity interval Ay = 3. 

and 158 AGeV. The study of charge fluctuations as a function of Ay was suggested in the 
original proposal J5|,|6[. The measured values of A$g vary between and —0.05. They are 
significantly larger than the values expected for QGP fluctuations (—0.5 < A$y < —0.15 



15|,|23|). The energy dependences of A$g for the largest rapidity interval {Ay = 3) and 



for the rapidity interval Ay = 1.2 are shown in Fig. [7| A weak decrease of A$q with 
increasing energy is suggested by the data. The numerical values of A$q for Ay = 1.2 
and Ay = 3 are given in Table ITTl Note that the fraction of the accepted tracks for a flxed 
rapidity interval Ay changes with collision energy, but this alone should not affect A$q 
provided the correlation length is smaller than the acceptance interval in y. 

VI. DISCUSSION 

The study of charge fluctuations in A+A collisions was motivated by the hypothesis 
that they may be sensitive to the creation of the Quark Gluon Plasma at the early stage 
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FIG. 5: The dependence of the measure <I>g on the fraction of accepted particles for central 
Pb+Pb collisions at 20-158 AGeV. Note that experimental points for a given energy are cor- 
related as the data used for a given rapidity interval are included in the broader intervals. 
The dashed line shows the dependence expected for the case when the only source of particle 
correlations is the global charge conservation, Eq. 3. 

of the collisions. To quantify the expected effect a simple QGP model was proposed in pj. 
In this model quarks and gluons are assumed to be in local equilibrium. Assuming entropy 
and net charge conservation during the evolution from the QGP to the final hadron state 
in each rapidity interval the number N of pions and their net charge is calculated. The 
number of charged pions is taken to be Nch = | ' ^ based on isospin symmetry. Using 
this model it was shown that the electric charge fluctuations are significantly smaller in 
the case of QGP creation than in the case of formation of confined matter at the early 
stage of the collisions (see Fig. |H)) p, |l5[- 
However, this model is not complete. A large fraction of pions originates from decays 



of resonances 



22J . This effect is expected to lead to a distortion of the charge fluctua- 



tions established after hadronization. To quantify this effect the model was extended as 
follows. From the total number of produced final state pions the entropy of the system 
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FIG. 6: The dependence of A$g on the width of the rapidity interval Ay for central Pb+Pb 
collisions at 20, 30, 40, 80 and 158 AGeV. Note that experimental points for a given energy are 
correlated as the data used for a given rapidity interval are included in the broader intervals. 

is calculated. This entropy is attributed to the early stage QGP, treated as an ideal 
gas of massless quarks and gluons. Bose-Einstein- and Fermi-Dirac-statistics are used to 
calculate equilibrium numbers of quarks and gluons. The rapidity distribution of these 
partons is centered at y = and is assumed to be of Gaussian shape with a = 0.8. For 
the calculations of charge fluctuations the rapidity interval —3 < y < 3 is divided into 
several (10 and 20) bins and in each bin the entropy and the net-charge of the contained 
partons is calculated. The resulting values of A$g at the QGP level are shown by the 
dashed line in Fig. 8. In the next step the QGP entropy is attributed to an ideal gas 
of p mesons. The numbers of p^, p^ and p° mesons in each bin are calculated assuming 
that I of all p mesons are neutral. Furthermore, all p mesons are assumed to decay into 
two pions. The rapidity distribution of the pions is divided into 20 bins and the A$q is 
calculated from the number of positively and negatively charged pions in each bin. The 
results of this model are shown in Fig. |Hlby the dotted curve. Note that by construction 
A$g = for the full acceptance. 
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FIG. 7: The energy dependence of A$q measured in central Pb+Pb collisions for a narrow 
rapidity interval Ay = 1.2 (left) and a broad rapidity interval Ay = 3 (right). 



TABLE II: The values of A$g for Ay = 3.00 and for Ay = 1.2 in central Pb+Pb collisions at 
20, 30, 40, 80 and 158 ^GeV. The first error is systematic, the second statistical. 



A$„ 



Ay = 3 
Ay =1.2 



20 AGeV 



-0.023±0.006±0.0001 
-0.023± 0.006± 0.0001 



30 ^GeV 



-0.028±0.0003±0.002 
-0.025±0.0002±0.016 



40 ^GeV 



-0.024±0.008±0.0005 
-0.016±0.008±0.0003 



A$. 



Ay = 3 
Ay = 1.2 



80 ^GeV 



-0.051±0.011±0.0002 
-0.040±0.011±0.0003 



160 ^GeV 



-0.036±0.013±0.0003 
-0.053±0.013±0.0004 



As expected, the decays of resonances strongly modify the initial QGP fluctuations. 
The value of A$q increases from values between —0.4 and —0.5 (the lower line in Fig. 
IH} to values close to zero (the upper line in Fig. IH)), the value characteristic for a gas 
of pions correlated by global charge conservation only. The model demonstrates that 
the distribution of charged particles in the detector acceptance is strongly distorted by 
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FIG. 8: The dependence of A^q on the fraction of accepted particles in central Pb+Pb collisions 
at 20-158 ^GeV. The prediction for the ideal QGP is indicated by the dashed curve (QGP), 
whereas the prediction for the QGP including hadronization and resonance decay is shown by 
the dotted curve (QGP+hadronization). 

the decay of intermediate resonance states. This may explain why the measurements do 
not show the suppression of the charge fluctuations naively expected in the case of QGP 
creation. 

The influence of resonance decays on charge fluctuations depends on the size of the 
rapidity interval Ay, in which fluctuations are calculated. If Ay is much bigger than 
the typical distance in rapidity of the daughter particles, the charge within the interval 
will not be changed by the decay and therefore the charge fluctuations should not be 
affected. On the other hand, if Ay is small, a large fraction of daughter particles will leave 
the interval and the initial net-charge will be significantly changed. The mean rapidity 
difference of two pions originating from decays of a p(770) meson is approximately 1 unit 
of rapidity. Therefore in order to minimize the decay effect the rapidity interval should 
be much larger than 1. However, this constraint is difficult to fulfill at SPS and lower 
energies because the rapidity distribution of all produced particles is not much broader 
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than 1. This explains an approximately constant value of A$q calculated within the 
QGP+hadronization model as seen in Fig. |H| A rapidity interval which is large enough 
to be unaffected by the influence of resonance decays would contain almost all particles 
produced in a collision. The net-charge in this interval would then reflect the number 
of participant protons and the fluctuations would be determined by fluctuations of the 
collision centrality and not the particle production mechanism. Thus at SPS energies the 
measured charge fluctuations are not sensitive to the initial QGP fluctuations. At very 
high energies (when the rapidity distribution of produced particles is significantly broader 
than 1) the charge fluctuations may be a valid signature of QGP creation. 

VII. SUMMARY 

Results on event-by-event charge fluctuations in central Pb+Pb collisions at 20, 30, 
40, 80 and 158 AGeV are presented in terms of the A$g measure. 

The measured A$g values are close to zero, as expected for a gas of pions correlated only 
by global charge fluctuations. This value is significantly higher than that predicted for 
the creation of a QGP and hadronization into pions with local conservation of entropy 
and net-charge. A model which incorporates intermediate resonances is described in this 
paper. Its results show that the decay of p mesons may easily increase the initial QGP 
charge fluctuations to A$q ^ thereby completely masking a possible QGP signal at 
SPS energies. 

The slightly negative value of A$g indicates correlations between positively and negatively 
charged particles beyond those from global charge conservation. The origin of these 
additional correlations may be the final state Coulomb interactions or quantum-statistical 
effects. 
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